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ABSTRACT 

In the framework of the Herschel-WISH key program, several ortho-H 2 and para-ILO emission lines, in the frequency range from 
500 to 1700 GHz, were observed with the HIFI instrument in two bow-shock regions (B2 and R) of the LI 157 cloud, which hosts 
what is considered to be the prototypical chemically-rich outflow. 

The primary aim is to analyse water emission lines as a diagnostic of the physical conditions in the blue (B2) and red-shifted (R) lobes 
to compare the excitation conditions. 

For this purpose, we ran the non-LTE RADEX model for a plane-parallel geometry to constrain the physical parameters (Ty„, Nh,o 
and n H2 ) of the water emission lines detected. 

A total of 5 ortho- and para-H2 16 plus one o-H2 ls O transitions were observed in B2 and R with a wide range of excitation energies 
(27 K<£' u <215 K). The H 2 spectra, observed in the two shocked regions, show that the H 2 profiles are markedly different in the 
two regions. In particular, at the bow-shock R, we observed broad (~30 km s -1 with respect to the ambient velocity) red-shifted wings 
where lines at different excitation peak at different red-shifted velocities. The B2 spectra are associated with a narrower velocity 
range (~6 km s _I ), peaking at the systemic velocity. The excitation analysis suggests, for B2, low values of column density N H2 o 
<5xl0 cm~", a density range of 10 5 <n H ,<10 7 cm , and warm temperatures (>300 K). The presence of the broad red-shifted wings 
and multiple peaks in the spectra of the R region, prompted the modelling of two components. High velocities are associated with 
relatively low temperatures (~ 100 K), Nh,o-5x10 12 -5x10 13 cut 2 and densities nH 2 -10 6 -10 8 cm 4 . Lower velocities are associated 
with higher excitation conditions with T kin >300 K, very dense gas (n H2 ~10 8 cur 3 ) and low column density (N H ,o<5xl0 13 cm -2 ). 
The overall analysis suggests that the emission in B2 comes from an extended (>15") region, whilst we cannot rule out the possibility 
that the emission in R arises from a smaller (>3") region. In this context, H 2 seems to be important in tracing different gas compo- 
nents with respect to other molecules, e.g. such as a classical jet tracer like SiO. We have compared a grid of C- and J-type shocks 
spanning different velocities (10 to 40 km s -1 ) and two pre-shock densities (2xl0 4 and 2xl0 5 cur 3 ), with the observed intensities. 
Although none of these models seem to be able to reproduce the absolute intensities of the water emissions observed, it appears that 
the occurrence of J-shocks, which can compress the gas to very high densities, cannot be ruled out in these environments. 

Key words. Stars: formation — Stars: low-mass — ISM: individual objects: LI 157 — ISM: molecules — ISM: jets and outflows 



1. Introduction 

Water controls the chemistry of many other species, whether in 
gas or solid phase and it is recognised to be a unique diagnos- 
tic of warm gas and energetic processes occurring in region of 
star formation (e.g., van Dishoeck et al. I201 11 1. H2O is a power- 
ful probe of the physical variations and the temporal evolution 
of outflowing material, which is processed by shocks produced 
by fast protostellar jets impacting against high-density mate- 
rial. Water, unlike most molecules, cannot routinely be observed 
with ground-based facilities. However, space instruments, such 
as SWAS and Odin, have provided the first detection of the 
ground state line ortho-FbO line at 557 GHz (e.g., Franklin et 
al. I20081 Bjerkeli et al. 120091 . The profile of this ground state 
line has been compared to that of CO line. This comparison is a 
direct testament to the H2O abundance being enhanced by many 
orders of magnitude in the shock (Kristensen et al. 2010). ISO 
was also capable of detecting a number of water transitions ex- 



cited in gas warmer than ~80 K (e.g., Liseau et al. 1996), but 
with spectrally unresolved data. 

The LI 157 bipolar outflow is the archetypical of the so- 
called chemically rich outflows (Bachiller & Perez Gutierrez 
[19971 hereafter BP97, Bachiller et al. l200Tb . At a distance of 
250 pc, and being favorably oriented in the plane of the sky, 
the LI 157 outflow is an ideal laboratory for observations of 
shocks chemistry. This outflow is known to be driven by a low- 
luminosity (~ 4 Lq) Class protostar and it is associated with 
several blue-shifted (B0, Bl, B2) and re d-shifted (R0, Rl, R, 
R2) lobes seen in CO (G ueth et al. 1199811, and in IR H 2 images 
(e.g. Neufeld et al. [19951 Nisini et al. 121)101 

The projected velocity of the northern (red-shifted) lobe is 
about 65% higher than that of the southern (blue-shifted) lobe 
(BP97) (see Fig. [TJ. However, the red-shifted lobe is more ex- 
tended than the blue-shifted one and the kinematical ages of 
the two lobes are found to be the same (15000 yr) indicating 
that both lobes were created simultaneously. Each outflow lobe 
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Fig. 1. PACS image of LI 157 of the integrated H2O emission 
at 1669 GHz (from Nisini et al. 120101 . Offsets are with respect 
to the L1157-mm source, at coordinates anooo = 20 h 39 m 06 s .2, 
£j2000 = +68° 02' 16'.'0. The contours of the SiO 3-2 integrated 
intensity is superimposed to the water map with a spatial reso- 
lution of 18" (Bachiller et al. I2001I) . The positions chosen for 
the HIFI observations (R and B2) are marked by two white cir- 
cles, which correspond to the largest and the smallest HIFI beam 
sizes. 

shows a clumpy structure and it seems that for each blue clump 
there is a corresponding red clump symmetrically located with 
respect to the central source. The symmetry in the individual 
clumps suggests that the outflow has undergone periods of en- 
hanced ejection and the curved shape suggests that there have 
been variations in the direction of the driving wind. In addi- 
tion, the morphology is clearly S-shaped indicating the presence 
of an underlying precessing jet. The brightest blue-shifted bow- 
shock, Bl, has been imaged at high-angular resolution revealing 
a clumpy arch-shape structure located at the top of an opened 
cavit y (Tafa lla & Bachiller [T995l Gueth et al. [19981 Benedettini 
et al. 120071 Codella et al. |2009l . L1157-B1 is well traced by 
species released by dust mantles such as H2CO, CH3OH, and 
NH3, complex molecules (Arce et al. 2008 ), and by SiO, which 
is the standard tracer of high-speed shocks (e.g., Gusdorf et al. 
120081 ) . In addition, LI 157-B1 has been observed with HIFI and 
PACS spectrometers onboard the Herschel Space Observatory 
(Codella et al. 2010, Lefloch et al. 2010): the preliminary results 
confirm the rich chemistry associated with the B 1 position, also 
showing bright H2O emission. 

The other bow-shocks have been less studied. As a conse- 
quence, it is still unclear how the physical and chemical proper- 
ties vary along the axis of the LI 157 outflow as well as between 
the blue- and red-shifted lobes. 



We present the analysis of various water emission lines ob- 
served with HIFI towards two shocked regions in the LI 157 out- 
flow, which is observed as part of the the key program WISH 
(Water In Star-forming-regions with HerscheQ, van Dishoeck et 
al. 120111) . The main aim of this paper is to constrain the physical 
parameters of the gas traced by water lines at the selected B2 and 
R positions of LI 157. We investigate the different excitation con- 
ditions of H2O in this shocked ambient gas in order to compare 
these conditions with other shock tracers. Since the abundance 
of water can be strongly affected by shocks, we may thus infer 
which type of shocks exist in the outflow regions traced by water 
emission. 



2. Observations 

Figure [T] presents the PACS map of the water emission at 1669 
GHz (from Nisini et al. 2010). The H2O map (green contour) 
is overlaid with contours of the emission from SiO (3-2) (ma- 
genta contour) (Bachiller et al. 120011) transitions. The water map 
exhibits several emission peaks corresponding to the positions 
of previously-known shocked knots, labelled as B0-B1-B2 for 
the south east blue-shifted lobe, and R0-R-R2 for the north west 
red-shifted lobe. For our water line observations we selected the 
shock positions B2 and R because they have different physical 
and chemical characteristics as shown in many previous works 
(e.g. Bachiller et al. [200Tb . Observations, from 500 to 1700 GHz, 
were carried out between May and November 2010 with the 
use of the HIFI instrument (de Graauw et al. 2010) in dual- 
beam switch mode with a nod of 3' using fast chopping. The 
HIFI receivers are double sideband with a sideband ratio close to 
unity. The data were processed with the ES A-supported package 
HIPE05.1 {Herschel Interactive Processing Environment, Ott et 
al. 2010 ) for baseline subtraction and then exported after level 2 
as FITS files in the CLASS9O/GILDAS0 format. Two polariza- 
tions, H and V, were measured simultaneously and then averaged 
together to improve the signal-to-noise. We checked individual 
exposures for bad spectra, summed exposures, fitted a low-order 
polynomial baseline, and integrated line intensities as appropri- 
ate. The main beam efficiency (77 m b) depends on frequency and 
is calculated as described in Roelfsema et al. (1201 11 ). ranging 
from to 0.75 to 0.71 in the 535-1670 GHz range. The absolute 
calibration uncertainty was estimated to be ~10%. At a velocity 
resolution of 1 km s _1 , the rms noise is 2-20 mK (T& scale) for 
frequencies less than 1113 GHz and 60 mK at 1669 GHz. 



3. Results 

A total of 17 emission lines were detected. Table Q] lists the ob- 
served H2O lines, associated with a wide range of excitation en- 
ergies (27 K<£u<215 K). In addition, a detection of o-H 2 18 
1 10 — lot in B2 will help us to constrain the optical depth. Also, 
in the observed spectral bands, additional molecular transitions 
have been detected in the B2 bow-shock, such as NH3 (lo-Oo), 
HCO + (6-5), CH 3 OH-E (ll U o~10i, 9 ), and C ls O (5-4). On the 
other hand, the I3 CO (10-9) line has been observed but not de- 
tected (see Table[2](. 



1 http://www.strw.leidenuniv.nl/WISH/. 

2 HIPE is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HIFI, PACS and SPIRE consortia. 

3 http://www.iram.fr/IRAMFR/GILDAS 
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Table 1. List of molecular species and transitions observed with HIFI in LI 157 B2 and R. Peak velocity, intensity (in scale), in- 
tegrated intensity (not corrected for beam efficiency, F m ), minimum and maximum velocities (V m! „, V max ), and linewidth (FWHM) 
are reported. 
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17 Frequencies and spectroscopic parameters have been extracted from the Jet Propulsion Laboratory molecular database (Pickett et al. 1998 i. 
h With respect to the ground H2O level lot for the oitho-H 2 level and 0oo for the para-H20 level, respectively. 
c At a spectral resolution of 1 km s _1 . 
''Gaussian fit when the observed profile is a gaussian-like. 
'Measurements taken at the dip of the emission line. 



3.1. Water profiles 

The B2 and R water spectra shown in Fig. [2] (here in units of 
antenna temperature, T&) were smoothed to a velocity resolution 
of 1 km s" 1 , to allow comparisons between all water transitions. 

The spectra measured in R are consistently fainter than those 
measured in B2 and the water transition at 1669 GHz appears 
to be barely detected at R. In the B2 spectra the profiles show 
a weak red emission and those lines with E u <80 K show ab- 
sorption at ambient velocity. The B2 spectra also show bright 
emission and peaks close to the systemic velocity while the bulk 
of the R emission is definitely red-shifted. The B2 spectra have 
a narrower velocity range than those observed at the R position. 
This could reflect a different shock velocity or simply be a con- 
sequence of a different inclination on the plane of sky. 

Interestingly, the 557 GHz profile observed at B2 is associ- 
ated with a much narrower wing (highest blue-shifted velocity, 
-9 km s ) than that observed in Bl (-25 km s~l), studied in the 
framework of the CHESS GT Key Project (see Fig. 3; Codella 
et al. 120101 Lefloch et al. 2010). This could be an indication 
of lower shock velocities at B2. On the other hand, B2 has a 
broader red wing than B 1 , supporting the idea that the profiles 
are strongly affected by geometrical effects. Further comparison 
of H2O profiles from observations obtained in B 1 and B2, would 
be instructive and will be performed when the full set of CHESS 
water data will be available. 

We note that the profiles in B2 narrow with increasing line 
excitation, while the line widths (FWZI) in R remain almost 
constant. In particular, the emission produced by transitions to 
the ground state level are those that appear broader (see Fig. [2]). 
These spectra are also those affected by the absorption dip and 
whit a higher S /N level. However, we cannot say if the narrow 
feature is a real trend. 

It is also notable that the spectra in R show multiple peaks 
that present a clear dichotomy with one peak at — 1-20 km s _I 
associated with transitions with £ u <60 K and a lower velocity 
peak at ~+10 km s associated with transitions with £ u >136 K. 
Surprisingly the high velocity emission, usually assumed to be 
emitted from regions close to the fast jet driving the shocks, is as- 
sociated with the low excitation emission lines. The dichotomy 
observed in the R bow shock has not been detected in any of 
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Fig. 3. Comparison of water emission lines at 557 GHz and NH3 
(lo-Oo) observed at LI 157-B1 (red profiles) and B2 (black pro- 
files). The vertical dashed line indicates the ambient LSR veloc- 
ity (+ 2.6 km s" 1 from C ls O emission; BP97) 

the molecules observed so far in LI 157. However, as already 
observed in B2, the difference in the line profiles at various 
excitation energies could be caused by a large amount of self- 
absorption in the low velocity range of those lines connecting 
with the ground state level. In fact, we actually observe a dip in 
the water emission line at 557 GHz at the systemic velocity. To 
explain such a striking difference (over 10-20 km s _1 ) between 
the line shapes, as that observed in R, the absorption would have 
to be caused by outflowing gas. It is evident that further multi- 
line analysis is required to test this assumption. It should also be 
noted that the same dichotomy has been observed in LI 148 at 
the R4 position by Santangelo et al. (1201 11 1, where they clearly 
show that this trend is real. 

Figure [4] shows the intensity ratio measured towards B2 and 
R of transitions selected to have different excitation but observed 
within similar HPBW, in order to minimise effects due to beam 
dilution. The ratios are plotted only at those velocities at which 
both lines have S /N >3. We note that, besides the obvious ef- 
fects due to the presence of the absorption dip, the 0-H2O (2 12- 
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Fig. 2. Molecular line profiles observed towards L1157-B2 and R: species and transitions are reported in the panels. The vertical 
dashed line indicates the ambient LSR velocity (+ 2.6 km s _1 from C ls O emission; BP97) 



loi)/p— H2O (ln-Ooo) water ratio seems to increase with veloc- 
ity, while in R it is clear that the P-H2O (2n-2o2)/p-H20 (ln- 
Ooo) water ratio decreases in both the low and high velocity com- 
ponents with respect to the peak at ~8 km s _1 . This effect, with 
an estimated line ratio error <30%, reflects the distinctive di- 
chotomy that was observed in the R bow-shock where water 
transitions with different energy excitations peak at two differ- 
ent velocities. 

In B2, the o-H 2 (1 io-loi)/o-H 2 ls O (lio-loi) ratio is used 
to estimate the H2O optical depth, T16, assuming that the two 
water isotopologues are tracing the same material with the same 
excitation temperature. To derive the o-H^O opacity, we need 
to use the 16 0/ I8 ratio for which we take the assumed local ISM 
value of 560 (Wilson & Rood |1994j l. The measured Ti6 value in 
the low velocity line wings is about 2, which is in agreement with 
our non-LTE excitation analysis. However, as it cannot be con- 
firmed that the o-H^O line and the o-H^O line are associated 
with the same excitation temperatures, any constraints placed on 



this value of opacity to infer the column density at B2 would be 
highly uncertain. 

3.2. Other line profiles 

Finally, the spectral set-up allowed us to observe also the NH3 
(lo-Oo) fundamental line: a tentative detection (3cr) of NH3 is 
found at R, while relatively bright emission is detected at B2 
as found in the nearby blue shock B 1 (see Fig [3] Codella et 
al. 1201 0). The ammonia emission, in the LI 157 outflow, will 
be analysed in a forthcoming paper. Table p2] lists the serendip- 
itous detections observed at the B2 bow-shock, HCO + (6-5), 
CH3OH-E (11 M0 -10i, 9 ), 13 CO (10-9) and C ls O (5-4), used to 
infer the water abundances. The B2 additional spectra are shown 
in Fig.|5](here in terms of antenna temperature, r a ). Spectra were 
smoothed to a velocity resolution of 1 km s , they all peak at 
the systemic velocity and are all associated with a relatively nar- 
row line width (< 6 km s -1 ) with respect to the water spectra. 
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Table 2. List of molecular species and transitions observed with HIFI in LI 157 B2 and R. Peak velocity, intensity (in scale), in- 
tegrated intensity (not corrected for beam efficiency, F m ), minimum and maximum velocities (V m! „, V max ), and linewidth (FWHM) 
are reported. 
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17 Frequencies and spectroscopic parameters have been extracted from the Jet Propulsion Laboratory molecular database (Pickett et al. 1998 l. 

b At a spectral resolution of 1 km s . 
'Gaussian fit when the observed profile is a gaussian-like. 




Fig. 4. Line temperature water ratios and water profiles in R and 
B2. Line temperature ratios, as function of velocity, are only 
shown where the S/N is larger than 3 in both emissions. 

4. Excitation Analysis 

We ran the RADEX non-LTE model (van der Tak et al. [2Q07T > 
with collisional rate coefficients by Faure et al. (2007) using the 
escape probability method for a plane parallel geometry in order 
to constrain the physical parameters (T^, Nh 2 o and «h z ) of wa- 
ter emission. An orfho/para=3 ratio was assumed, equal to the 
high temperature equilibrium value. To place some constraints 
on the unknown size of the emitting region, we assumed three 
different sizes (3", 15" and 30"). RADEX does not take into ac- 
count the near-IR pumping of the H2O lines where the radiation 
field is relatively strong and has some impact on the excitation 
conditions of the emitting gas. However, since both B2 and R 
position are relatively far from the central source, we expect no 
continuum to affect our results. 

4.1. L1157B2 

At B2, given the presence of the absorption dip and the rel- 
atively narrow velocity range, we model the physical param- 
eters of water emission by integrating the H2O intensity over 
the whole spectral profile (i.e. neglecting the absorption dip). 
We will consider the measurements of water transitions with an 



13 C0 (10-9) 




-20 -15 -10 -5 5 10 15 20 25 30 
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Fig. 5. Additional emission lines (beside the H2O and NH3 pro- 
files reported in Fig. [2]) observed at LI 157-B2. 

absorption dip as upper limits. In Fig. [6] we overplot non-LTE 
RADEX model predictions against integrated flux water ratios. 
For this analysis, we conservatively assumed an uncertainty of 
20%. Coloured triangles with uncertainties are the observed val- 
ues corrected for beam filling and assuming three different sizes 
as mentioned in the insert of the plots. The top panel of Fig. [6] 
shows a plot of the 0-H2O (3i2-3o3)/p-H20 (ln-Ooo) ratio 
versus 0-H2O (3i2-3cb), the middle panel shows a plot of the 
0-H2O (2i 2 -loi)/o-H 2 (lio-loi) ratio versus o-H 2 (2i 2 - 
-loi), while the bottom panel shows a plot of the 0-H2O (3i2- 
-3 03 )/o-H 2 O (2 12 -l i) versus o-H 2 (3 12 -3 3)/o-H 2 O (li - 
-loi) ratios. 

Some general conclusions can be drawn from the inspection 
of Fig. [6] The observations are consistent with model predictions 
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Table 3. Physical conditions traced by water lines. 
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only assuming an emission size >3" (black triangle). In particu- 
lar, a size closer to 15"— 30" in agreement with the LI 157 map 
of the o-H 2 (2i 2 — lot) line published by Nisini et al. (120101 1 
is suggested (see bottom panel). The observed water ratio ap- 
pears to be well constrained from low values of column density 
Nh,o <5x10 13 crrT 2 (see the top panel), while a density range of 
10 5 <nH 2 <10 7 ctrT 3 , and temperatures Tki n >300 K are inferred 
from all panels. However, because of the presence of the absorp- 
tion dip in both water transitions plotted in the middle panel, 
the inferred physical conditions can be only considered as upper 
limit. In fact, the absorption dip could affect our analysis, which 
would shift observations towards lower excitation conditions. 

With regards to the absorption dip observed in the water 
ground state transition, if we assume that the absorption is due 
to foreground gas unrelated to LI 157-B2, we estimate, from the 
standard radiative transfer equation, a 8 K<r ex <13 K adopting a 
0.1<Ti6<10. This low value of T ex usually implies low densities 
or low kinetic temperatures. No other constraints can be inferred 
due to the uncertainties in the measurements and too many free 
parameters to fit simultaneously with the RADEX code. 

4.2. L1157R 

In the water profiles of the bow-shock R, Fig. [2]shows an evident 
dichotomy at different excitations where two emission peaks at 
different velocities clearly displaying different excitation condi- 
tions. Thus, we modelled the emission by splitting the lines into 
two components (hereafter called LV, at — hlO km s _1 , and HV, 
at ~+20 km s" 1 ). In Fig. [7] we overplotted non-LTE RADEX 
model predictions against the observed LV and HV T m \, water 
ratios. The top panel of Fig. [7] shows a plot of the 0-H2O (3i2- 
-303V0-H2O (1 10-I01) ratio versus 0-H2O (3i2-3o3), while the 
bottom panel of Fig [7] shows a plot of the P-H2O ( 1 1 1 -0oo)/o- 
-H 2 (3i2 -3 3) ratio versus p-H 2 (ln-Ooo)- 

A first conclusion that can be drawn is that the LV compo- 
nent appears to be associated with higher excitation conditions 
than the HV component and we can exclude a point-like source 
(i.e. the 3" case) based on the inconsistency of column density, 
which can be derived from the two panels. For the low excitation 
HV, we have standard kinetic temperatures >100 K and a de- 
generacy between column density and hydrogen density: either 
Nh 2 o~5x10 13 cm" 2 and n H2 ~10 6 cirT 3 or N H2 o~5xl0 12 cirT 2 
and ne,~10 7 cm 3 . On the other hand, for the high excitation LV 
emission the observed water ratio traces temperatures Tk m >300 
K, as well as surprisingly dense gas with nH,~10 8 cirT 3 and low 
column density, Nh,o<5x10 13 cirr 2 . The tentative detection of 
NH3 (lo-Oo) emission at these velocities seems to support a very 
high density solution. 

5. On the origin of water emission in L1 1 57 

In summary, Table [3] outlines the physical conditions traced by 
water that we constrain from our analysis. The water at B2 traces 
lower density environments than that at R-LV and R-HV. In addi- 
tion, the main difference between R-LV and R-HV is that H2O at 
R-HV traces lower temperatures and densities with respect to R- 



LV. The physical parameters derived in B2 and R-HV are in rea- 
sonable agreement with those measured by Nisini et al. (2007), 
using multi transitions observations of a classical jet tracer like 
SiO. On the other hand, the density inferred for R-LV (nn,~10 8 
cm -3 ) is higher with respect to that derived from SiO emission 
by two orders of magnitude. Actually, the SiO profiles as ob- 
served at the R position (see Fig. |4]l cast serious doubts on the 
theory that emission from these species has a common physi- 
cal origin. In addition, our results strongly support the conclu- 
sion of Nisini et al. (2010) that water emission seems to follow 
quite closely the distribution of H2 emission. In fact, the general 
physical conditions inferred in the R bow shock (Tk m >300 K, 
nH,~10 6 -10 7 cm 4 ) are consistent with those derived for the H2 
in Nisini et al. d20T0l l. 

5.1. Water abundances 

We derive, at the R bow shock position, the H2O abundance, 
with respect to H2, using the 7V(H2) estimated by Nisini et al. 
(2007) and in B2 using the C ls O (5-4) line emission detected in 
the present survey (see Fig.|5]and Tab.|2|. The velocity averaged 
H 2 abundance for R is -lO^-lO" 7 . 

The H2O abundance for B2 is estimated to be ~10~ 6 assuming 
[C 18 O]/[H 2 ]=2xl0- 7 (Wilson & Rood fT994l l in the temperature 
range from 300 to 500 K (see Table |3). This result is consistent 
with the water abundance found at B 1 at low velocities (Lefloch 
et al. 1201 Oi l. On the other hand, the X(H 2 0) measured at the 
highest velocities in Bl is 2 orders of magnitude greater than 
in B2 (~10~ 4 ). These differences could provide evidence for an 
older shock in B2 compared to the one in B 1 . 
The results found in R are comparable with those obtained by 
Santangelo et al. (1201 11 1 in the L1448 outflow. These low wa- 
ter abundances could support the possibility of having a J-shock 
instead of a C-shock where the H2O abundance is expected to 
be ~10~ 4 relative to H2. Alternatively, the low water abundances 
measured could also be evidence of either an old shock where 
the water in the gas phase has had time to be depleted on grains 
or of UV dissociation of H2O (see Bergin et al.;1998). 

5.2. Comparison with shocks models 

Theoretical studies indicated that, under typical interstellar con- 
ditions, shocks faster than ~20 km s _1 are efficient enough to 
free much of the water ice frozen on grain surfaces (Draine et al. 
1983), while shocks faster than ~15 km s _1 , through gas-phase 
chemical reactions, produce large quantities of water (Kaufman 
& Neufeld I19961 I. After the shocked gas has cooled enhanced 
water abundances of ~10~ 4 relative to H2 can endure for as long 
as 10 5 yr (Bergin et al. 1998). Sub-millimeter and far-infrared 
water emission lines are very efficient coolers, therefore large 
interstellar water abundances produced by the shock processing 
are relevant for the thermal evolution of the gas (e.g., Neufeld et 
al. 119951 1. Observations of a number of excited molecular tran- 
sitions in LI 157 clearly indicate that shocks are present in both 
lobes. The present dataset deserves a more detailed comparison 
with shock models. In this present paper we only discuss our 
data using models provided by literature. 

We have compared a grid of C- and J-type shocks spanning 
different velocities (10 to 40 km s~') and two pre-shock densi- 
ties (2xl0 4 and 2xl0 5 crrr 3 ) provided by Flower & Pineau des 
Forets (120101 1. with the observed intensities. Figure [8] shows the 
observed line ratios with respect to the high energy (£' ll =215 K) 
1097 GHz line for different sizes of B2, R-LV and R-HV against 




Fig. 6. Analysis of the water line emissions in L1157 B2. 0-H2O (3i2-3o3)/p-H20 (In— Ooo) ratio versus 0-H2O (3 12— 3o3) (in 
logarithmic scale, top panel), 0-H2O (212-I01V0-H2O (lio-loi) ratio versus 0-H2O (2i2-loi) (middle panel), and 0-H2O (3i2- 
-3o3)/o— H2O (1 10— lot) versus 0-H2O (3i2-3o3)/o-H20 (2i2-loi) ratios (bottom panel), for non-LTE (RADEX) plane parallel 
models at the labelled temperatures. A linewidth of 5.5 km s has been assumed, according to the observed spectra. Each coloured 
curve corresponds to the labelled Nh 2 o (see insert at the top-left corner). H2 density increases from left to right as the labelled 
values. Triangles with error bars indicate the size of the emitting region (30", 15", 3"). 



the predicted shock line ratios. From our analysis, none of these 
models seem to be able to properly reproduce the absolute inten- 
sities of the water emissions observed. 



Particularly in the case of the R-HV component, a degen- 
eracy is noticeable because, for an emitting region of size 30", 
both a speed of 10 km s with a pre-shock density of 2xl0 5 
cm -3 and a speed of 30 km s _1 with a pre-shock density of 2x1 4 
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Fig. 7. Analysis of the water line emission in L1157 R. 0-H2O (3r2-3o3)/o— H2O (1 10— lot) ratio versus 0-H2O (3i2-3o3) (in 
logarithmic scale, top panel), and P-H2O (1 n-Ooo)/o-H20 (3i2-3(b) ratio versus P-H2O (ln-Ooo) (bottom panel), for non-LTE 
(RADEX) plane parallel models at the labelled temperatures. A linewidth of 10 km s has been assumed, consistent with the 
observed spectra. Each coloured curve corresponds to the labelled A^h,o (see insert). H2 density increases from left to right as the 
labelled values. Triangles with error bars indicate the size of the emitting region (30", 15", 3"). 



cm can be inferred for both types of shock. This is consistent 
with the analysis made in the water diagnostic plots (Fig.|7]i. The 
same sort of degeneracy can be drawn for the R-LV component, 
where we do not find a satisfactory solution that matches the 
full sample of water transitions. The most plausible explanation 
is given by, in agreement with Fig. [7] a size of 3" for a C-type 
shock with a speed of 30 km s _1 and pre-shock density of 2xl0 4 
cm -3 , and for a J-type shock with a speed of 20 km s" 1 and pre- 
shock density of 2xl0 5 cm 4 . 

Regarding the B2 position, a lower speed of J-shocks around 
10 km s is inferred. Two different solutions are found: (i) a 
pre-shock density of 2xl0 4 and 3" of size, or (ii) a pre-shock 
density of 2x 10 5 and 15" of size. According to our previous wa- 
ter analysis an emitting size of 15" looks more plausible. These 
preliminary results call for detailed modelling for a more com- 
plete analysis. Note, for instance, that UV H2O dissociation is 
not taken into account in the grid of shock models used. 

C-shocks have been recently invoked by Flower & Pineau 
des Forets d20T0j l and Viti et al. (f20TTb modelling LI 157 Bl. On 
the other hand, recent results obtained with PACS in the course 
of the CHESS spectral survey key-program (Benedettini et al. in 
prep.) required J-shocks to explain the CO excitation conditions 



in LI 157 Bl. This scenario could also provide a possible expla- 
nation of our analysis in B2 and R where both C- and J -type 
shocks seem to be required. 

6. Conclusions 

We have presented a water emission observed in the 500-1700 
GHz band towards the two bow-shocks B2 and R, of the LI 157 
proto-stellar outflow. The main conclusions are the following: 

1. The comparison between H2O and SiO profiles, and their 
physical characteristics casts serious doubt on the assump- 
tion that emission from these species has a similar physical 
origin. 

2. We derive H 2 abundances for R -lO^-lO 7 and for B2 
~10~ 6 . This result is consistent with the water abundance 
found at B 1 at low velocities. On the other hand, the X(H20) 
measured at the highest velocities in Bl is 2 orders of mag- 
nitude greater than in B2 (~10~ 4 ). These differences could 
provide evidence for an older shock in B2 compared to that 
inBl. 

3. The emerging scenario highlights the importance of J- 
shocks, expected to be associated with a thin layer and very 
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Fig. 8. Comparison between the observed HIFI H2O line ratios of the two velocity components in R and B2 for different sizes (black 
triangles) and those predicted by C-type (red triangle) and J-type (cyan triangles) shock models from Flower & Pineau des Forets 
(l2010l l. 



densely compressed material in these environments. The re- 
sults found in R are comparable with those obtained by 
Santangelo et al. (submitted) in the L1448 outflow. These 
low water abundances could support the possibility of hav- 
ing a J-shock instead of a C-shock. 
4. Interestingly, the highest excitation conditions are observed 
at low velocities. However, if we assume that high excitation 
is tracing portions of gas near the shock, we could solve this 
apparent contradiction by assuming that we are observing a 
very collimated region located along the plane of the sky. In 
this case the fast collimated gas should be reprojected and 
thus have lower radial velocities. On the other hand, a less 
excited region could be more extended, which would result 
in a wider velocity range due to the geometry. 
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